Cumulative compressibility effects on population dynamics in turbulent flows 
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Bacteria and plankton populations living in oceans and lakes reproduce and die under the in- 
fluence of turbulent currents. Turbulent transport interacts in a complex way with the dynamics 
of populations because the typical reproduction time of microorganism is within the inertial range 
of turbulent time scales. In the present manuscript we quantitatively investigate the effect of flow 
compressibility on the dynamics of populations. While a small compressibility can be induced by 
several physical mechanisms, like density mismatch or the finite size of microorganisms with respect 
to the fluid turbulence, its effect on the the carrying capacity of the ecosystem can be dramatic. 
We report, for the first time, how a small compressibility can produce a sizeable reduction in the 
carrying capacity, due to an integrated effect made possible by the long replication times of the 
organisms with respect to turbulent time scales. A full statistical quantification of the fluctuations 
of population concentration field leads to data collapse over a broad range in parameter space. 
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Bacterial colonics living on a nutricnt-rich hard agar 
dish are commonly described in terms of the Fisher- 
Kolmogorov-Petrovsky-Piscounov (FKPP) equation of 
population dynamics [l[ . On a rigid substrate the FKPP 
equation includes only diffusion and growth terms. In 
a fluid environment the flow can play an important role 
by means of the advective transport of microorganisms, 
both at low [ij, as well as at high Reynolds numbers. 
In particular, in oceans and lakes microorganisms have 
found ways to thrive and prosper in high Reynolds num- 
ber fluid environments. 

As discussed in 0, 01, the advection by means of 
a compressible turbulent flow can lead to highly non- 
trivial dynamics. Two striking effects emerge: for small 
enough growth rate, /i, the population concentration 
field, c{x,t), is strongly localized near the transient but 
long-living sinks of the turbulent flows; in the same limit, 
the average space-time concentration of the population 
(carrying capacity) becomes much smaller than its max- 
imum value in absence of flow. Both effects have impor- 
tant biological implications. Recently it was shown that 
this phenomenology holds true in the case of more realis- 
tic two dimensional compressible turbulent velocity fields 



In this manuscript we present new results that support 
the key role played by even a slight level of compressibil- 
ity. We have in mind a simplified model for photosyn- 
thetic microorganisms that actively control their buoy- 
ancy to maintain themselves at a fixed depth below the 
surface of a turbulent fluid. Small density mismatches 
or inertial effects due to the finite size of the microor- 



mnisms can also lead to an effectively compressible flow 



The main result of our investigation is that the overall 
carrying capacity can be described as a universal function 
of the non-dimensional growth rate and compressibility. 
Our results imply that even a small compressibility leads 
to a substantial reduction in the overall carrying capacity. 

Bacterial populations in nutrient rich environments 
can be modeled in terms of the following Fisher- 
Kolmogorov-Petrovsky-Piscunov (FKPP) equation: 

^ ~ (uc) = DV^c + ficil- c). (1) 



The above equation describes the evolution of the mi- 
croorganisms concentration field, c(x, i), under the effect 
of advective transport by means of a (turbulent) velocity 
field, u{x,t), and in presence of spatial diffusivity with 
a diffusion constant, D, and a replication rate, fj,. The 
parameters D and /i encode important biological infor- 
mation which may depend on the amount of nutrients, 
temperature, organismic motility and many other param- 
eters of the ecosystem. In Eq. ([T]) the averaged popu- 
lation density in absence of flow which can be sustained 
by the ecosystem has been rescaled to 1; hence our con- 
centration field c is dimensionless. As an example of "life 
at high Reynolds number" one could consider the FKPP 
equation [Eq. ([T])] to represent the density of the ma- 
rine cyanobacteria Synechococcus under conditions of 
abundant nutrients (so that /i ~ constant) . 

In this study we build a simple model of microorgan- 
isms living in a localized layer at some depth e.g. un- 
der the ocean surface. For this reason we will consider 



2 



e 






E e 















0.0 




3.17 0.61 


5.80 


■ 10" 


2 


1.0 




15 


1.45 ■ 10" 


2 


3.14 0.55 


6.10 


■ 10" 


2 


6.8- 10" 


1 


30 


6.50- 10" 


2 


2.76 0.47 6.60 


■ 10" 


2 


3.2- 10" 


1 


45 


1.62 • 10' 


1 


2.24 0.36 


7.60 


■ 10" 


2 


1.6- 10" 


1 


60 


3.90 • 10" 


1 


1.68 0.25 


8.90 


■ 10" 


2 


4.1 • 10" 


2 


75 


7.50 ■ 10" 


1 


1.38 0.15 


1.16 


■ 10" 


1 


8.4- 10" 


3 


90 


1.0 




1.25 0.12 


1.31 


• 10" 


1 


1.6- 10" 


3 



TABLE I. For different values of 6, the compressibility 
the energy E = J u^dxdy, the energy dissipation rate 
e = Yi J {'Vu)^dxdy, the Kolmogorov time scale r,, used in 
our simulations. Different values of the compressibility are 
obtained by projecting out the compressible and incompress- 
ible component of our 2d compressible field for k = 0.16 as 
discussed in the text. The energy and dissipation are rounded 
off to two decimal places. 

a two-dimensional planar surface, taken out of a fully 
three dimensional turbulent flow. A projected two di- 
mensional slice from a fully 3d velocity field is compress- 
ible. The dimensionless numbers characterizing the evo- 
lution of the scalar field, c{x, t), are the Schmidt number 
Sc = v/D and the dimensionless growth rate /zr,,. Here 
T,, = [yjef-^'^ is the Kolmogorov dissipative time scale 
and £ is the energy dissipation rate of the fluid. A par- 
ticularly interesting regime arises when the growth time 
Tg = yU^^ lies in the middle of the inertial range. Al- 
though many factors can affect the estimates, this is the 
case for oceanic cyanobacteria and phytoplankton living 
in an turbulent ocean. Indeed, oceanic turbulence has ed- 
dies with turnover times ranging from minutes to months 
while typical microorganism doubling times range from 
half an hour to half a day 0, II] . 

Details of the numerical simulations can be found in 
j^]. The velocity time series of the two dimensional 
slice were filtered in order to produce time histories 
with different degrees of compressibility. By projecting 
the velocity field into an incompressible, m\ and com- 
pressible part, m'^, we could define a family of velocity 
fields u = a/2[m'= sin(6') + cos(6')] with < 6* < tt/2. 
The dimensionless compressibility of each new flow is 
K = (V • ii)^/(Vw)2. For 6* = the velocity field is 
purely incompressible, k = 0, whereas for = 7r/2 it is a 
potential flow, k = 1. The case ~ ii/A corresponds to 
the flow studied in Q ■ In Table U we report the values 
of the compressibility used. The carrying capacity of the 
ecosystem is defined as: 

Z={^j c(x,t)d2x), (2) 

where the integral runs over a square domain of size 
with periodic boundary conditions and () indicates tem- 
poral averaging. In Fig. [T] we report Z as a function 
of the dimensionless growth rate, /zr^, for different val- 
ues of the flow compressibility (we recall that in absence 
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FIG. 1. Time-averaged carrying capacity Z versus the non- 
dimensional growth rate /xt^ for different values of compress- 
ibility hi. Region with /jr,, < TnjrL and [lT^ > 1 has been 
shaded for visual aid. Here, tl = L/y/E is the box size time 
scale with L = 27r 0| . 

of flow Z = 1). The shaded region on the right repre- 
sents /iT,, > 1, while the left shaded region reflects the 
condition /xr,; < r^/r^, where = L/\/E is the large 
eddy turnover time. For large growth rates the curves 
saturate towards unity, although for smaller compress- 
ibility the drop in carrying capacity shifts to smaller and 
smaller growth rates, jUt,,. Note that even at very small 
values of the compressibility (e.g. k = 0.0145) there is 
a significant reduction of the carrying capacity, ^ 20%, 
for small growth rates. This regime is particularly rele- 
vant to marine biology where compressibility is small but 
the orgamisms have a very slow reproduction rate when 
compared with turbulent time scales. As a result, despite 
small compressibility, one should still expect important 
effects on the global carrying capacity. For large values 
of the growth rate, i.e. Tg <^ t^i, the time-scales of tur- 
bulence are too slow to have any effect on the evolution 
of the concentration, and Z — >■ 1. 

In the limit Tg ^ 0, the concentration field tends to 
become uniform with the leading correction coming from 
the local compressibility. After a series expansion one 
obtaines (Z)^ w 1 - iT^/L^){JiV ■ u)^dj<.) + 0{t^). The 
limiting values for small growth rates are complicated to 
access numerically, especially for k < 0.16, as it takes a 
longer and longer time for the carrying capacity to reach 
a steady state at decreasing fjL = t~^ . To obtain reliable 
estimates in this limit we can proceed as suggested in [3| 
where it was shown using a multifractal analysis that as 
/i — >■ the statistics of c closely resembles the statistics 
of the corresponding passive scalar probability density 
P{x,t), which satisfies Eq. ([1]) with /i = 0. The limiting 
values of Z^^^q for different values of compressibility is 
given in Table HI 

The spatial behavior of c{x,t) is characterized by 
strong fiuctuations with non trivial correlations in both 
space and time. In Q it was shown that the multifractal 
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FIG. 2. Log-log plot of the anomalous scaling exponent 02 
as a function for the compressibility n in the case = and 
/i = 2. Note that the two scaling laws have the same slope 
with different offset. Inset: Scaling behavior of {Cfi{r)) for 
fj, = Q and k = 0.16. 

analysis describes the statistical properties of the con- 
centration field. Here wc extend the analysis to different 
values of the compressibility, k, and of the growth rate, 
ji. We focus on the scaling of: 

^mW = ^(/ c{x,tYd?x) (3) 

r \JB(r) I 

where (...) now stands for a time and ensemble aver- 
age and B{r) is a disk of radius r. Within the inertial 
range, we find that C^{r) ^ r^^i with Ooa non linear 
function of q constrained to og = ai = ]^, [l3| (in the 
inset of Figure [5] we show the scaling of (r) versus r 
for a particular value of /i = and k = 0.16; plots for 
other values look similar). It is important to remember 
that each k > can be associated with a compressibility 
length (and time) scale, Ik = u„ns/((V • m)^)^/^ where 
Urms = is the root-mean-square velocity. For very 
small compressibilities, this scale will be comparable or 
larger than the integral scale of turbulence and the sys- 
tem can effectively be considered as incompressible. In 
our system this happened roughly for values of k < 0.065 
and indeed for these compressibility values no scaling can 
be detected, making it impossible to compute a^. Here- 
after we limit our analysis to the region k > 0.065. 

We expect the scaling exponent a2 to be a function of 
/i, r,, and k. Using dimensional analysis, we expect that 
02 depends on fi and r,, through the dimensionless com- 
bination ^Trj and that the time scale associated to the 
velocity gradient in the turbulent flow with compressibil- 
ity k, should be given by l/((V-u)-^)-^/^ = T,,/-y/K. Hence, 
we conjecture that a2 should depend only on the dimen- 
sionless combination ^JUj [iTr^. Since for /i = we find 
that a2 has a finite value for any value of compressibility 
K, we further conjecture that a-i cx y^/(/^Tj, -|- jS) where /3 
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FIG. 3. Logarithmic plot of the carrying capacity Z versus 
the anomalous exponent a-i,(K) for ^ = and fi — 2. The 
black line represent a best fit with slope —3.9 ± 0.1. We fit 
the data separately through jj, — and /j, = 2 datasets. The 
mean of the fits is the slope and the error is the maximum 
deviation of the fit from the two slopes. Inset shows the linear 
relationship between the localization length and the estimated 
cut-off length for Sc = 0.1,0.5,1, and 10. The black line 
represents the best fit and has a slope ~ 0.6. 

is a constant. We expect that /3 « t^/tl ~ l/Re\, where 
tl is the large eddy turnover time discussed above. 

In Fig. [21 we show the scaling exponent 02 for /.i = 
and /i = 2 and for different compressibilities n. The most 
striking feature of Fig. [2] is the well defined scaling law 
between 02 and n and that the scaling exponent does not 
change by changing /^r^. Therefore we can write: 

a2(K,/^T,,) = a(pTn)K.'' (4) 

where a(/XT^) = a/{^Trj + (5). Using the data in Fig. [2] 
we obtain 7 = 0.62, reasonably close to the prediction 
7 = 0.5 given by dimensional analysis. Note that the 
function a(/iT^) in Eq. Q is a decreasing function of /ir^, 
and that lim^jT-^^o 02 (//r^) = 02(0) > 0. 

In Ref. [3l it was shown that for ^ ~ the carrying 
capacity is related to 02 by Z w {^q/LY^ where = 
{P'^)/{{VPf) and P{x,t) is the solution of Eq. ^ with 
/i = 0. Here we generalize the results of Ref.[3| for /i 7^ 
and compressibilities n (see Table. |l|. In Fig. [3] we plot 
log(Z) versus a2(/iT^,K) for /i = and /.i = 2, which 
supports the scaling ansatz 

Z{fI,Tr„K) = i^j (5) 

which defines the length scale ^est- In order to gain a 
deeper physical insight into the meaning of the cutoff 
scale £,est we define a scale based on the gradient of the 
concentration field which is a generalization of our pre- 
vious definition of 
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FIG. 4. Semi-logarithmic plot of Z as a function of {fJ-Tn + 
P)/^? where 7 = 0.62 is given by the fit of Fig.[2]and /3 = 0.1. 
All the data of Fig. [1] has been collapsed on a single curve. 



In the inset of Fig. |3] we show the behavior of the two 
length scales ^est and by varying the Schmidt number 
0.1 < Sc < 10. Evidently the two definitions are con- 
sistent within a numerical prefactor, suggesting that ^est 
is proportional to the cutoff scale of the gradients of the 
concentration field. 

We can now test for a universal scaling behavior of 
the carrying capacity for varying compressibility k and 
non-dimensional growth rate ^r,, using Eq. [5l 

Upon defining the quantity X = log(a(/ir^)K''') the 
validity of Eq. ([5]) can be checked by replotting the 
data of Fig. [1] for Z as a function of X. Since X = 
log(Q;) + \og[K< /{^1T^ + /3)] and 7 = 0.62 from Fig. 
data collapse of the results shown in Fig. [T] should hold 
for some particular value of /3 whereas the exact value 
of a represents just a horizontal shift in the whole data 
set. In Fig. |4] we show that for /? = 0.1 all the data 
collapse on a well defined curve which represents our 
universal function. This estimate /3 = 0.1 obtained for 
the data collapse is very close to the dimensional esti- 
mate /3 ~ Tri/TL ~ 0.02 — 0.03 (shaded region on the 
left in Fig. [ij. This result is important since it shows 
the deep and highly non-trivial connection between a 
bulk property of the system, namely Z , and the inter- 
mittency parameters of the FKPP equation in compress- 
ible flows. Moreover, Fig. |4] represents a prediction of 
the average carrying capacity covering an entire param- 
eter space spanned by the three basic quantities of the 
system, namely ^, k and t^. Our numerical results can 
thus be extrapolated to new regimes, to investigate the 
importance of weak compressibility due, for instance, to 
inertial effects or buoyancy forces in populations subject 
to oceanic turbulence. 

In summary, we have used the FKPP equation to nu- 
merically study the population dynamics in a compress- 



ible turbulent velocity field. As a simplified model rele- 
vant for marine biology we have considered microorgan- 
isms (bacteria or plankton) confined to a two-dimensional 
plane experiencing the effect of a three dimensional ve- 
locity field. We have investigated in detail the effect of 
compressibility in a velocity field over a wide range of k's. 
We found that even a small compressibility can signifi- 
cantly reduce a global quantity like the average carrying 
capacity, due to the slow reproduction rate of the or- 
ganisms. We expect that in oceans or lake this may be a 
common situation. We further quantified in terms of spa- 
tial intermittency exponents the statistical properties of 
the concentration field. Our study clearly suggests that 
it is quantitatively wrong to neglect even small degrees 
of effective compressibility. This compressibility, even in 
absence of organisms with active buoyancy control, can 
be induced by density mismatches or by the finite size of 
the organisms. Experimental tests of our findings would 
be of primary importance. 
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